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As a step in the direction of "computer vision," several programs 
have been written which transform the output of a vidissector into some 
mathematical descriptions of the boundaries enclosing the objects in the 
field of view. Most of the discussion concerns the techniques used to 
transform a sequence of points, presumably representing a curve in the 
two-dimensional plane of view, into the best-fit conic-curve segment, or 
best-fit straight line. The resultant output of this stage is a list of 
such segments, one list for each boundary found. 
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As a step in the direction of "computer vision,'' several programs 
have been written which, applied successively, transform the point-view 
from the vidissector into a mathematical description of the boundaries 
of the objects in the field of view. An executive program will co-ordinate 
the actions of the various programs described below in order to find each 
and every object in the view. For example, several attempts may have to 
be made in order to be sure that a single object has been seen, using such 
clues as color contiguity, graininess, or texture, light intensity, etc. 

The programs REGI0NS1 and BNDS0RT1 examine the points in the current 
subview of the vidissector and determine which ones belong to the object 
or region being sought « The executive determines which section of the 
vidissector to look at 5 the density of points to be interviewed, and the 
predicates that define a region. The end result is an ordered list of 
points consituting the boundary of the region. (Each point is a list of 
two LISP numbers, which are the x- and y -coordinates on the vidissector 
screen of the geometric point represented.) Further information on these 
routines may be obtained from Gerald Sussman. 

POLYSEG takes the output of BNDS0RT1 and a straight-line approxima- 
tion to the boundary of the figure, using techniques adapted from the work 
of Roberts (see Lincoln Laboratory Technical Report No. 315, n Machine Per- 
ception of Three-Dimensional Solids") . The purpose of this routine is 
twofold: (1) to smooth out minor pieces of noise in the data, and (2) to 
reduce the total number of points representing the boundary. Hopefully 
the error tolerance will insure that a sufficient number of short straight 
lines are used to approximate a curved line segment on the boundary. 
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Perhaps more work could be put into this routine in the form of a better 
heuristic determination of corners as they appear in a boundary. POLYSEG 
returns a list of n+1 points for an n-sided polygonal approximation to 
the boundary of the region. The first point of the list is duplicated 
as the last point also, thus explaining the additional point. 

ANGLES takes the output of POLYSEG and returns a list of four items 

(1) Augmented polygon list [see below] , 

(2) Reasonable approximation to the diameter of the region described 
by the input polygon list, 

(3) "CLCK" or "CNTK" indicating clockwise or counterclockwise ordering 
of the boundary points, 

(4) Area of the approximating polygon (positive, . floating point). 
The original polygon list, i.e. the input list, has only the x-and y-coor- 
dinates for each vertex; the augmented list, (1) above, provides also the 
inside angle at each vertex and the distance to the next-in-order vertex. 
Thus each member (x.y.) of the input becomes (x.y.a.d.) in the augmented 
output list. 

SLPCS takes the output of ANGLES (providing the points are ordered 
clockwise) and attempts to segment the boundary into true straight lines 
or curved pieces. By "piece" we mean a segment of the boundary which, 
according to the following criteria, is a straight line, or a straight- 
line approximation to some curve without inflections. The output is then 
a listing of the pieces of the boundary in order, and the listing of the 
points of each peice is ordered the same as in the input. In addition to 
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the points of the relevant boundary segment, the first member of each 
piece (which is otherwise a list of points) is a fixed point integer 0, 
1, or 2 corresponding to "straight line," "curved out," and "curved in" 
respectively. Example: 

1 

/ \ 2 








» 




Notice that each piece is complete in itself so that its endpoints are 
the same points as will be found on one end of either adjoining piece. 
Several heuristics are applied to decide where one piece should end and 
another begin. Clearly an inflection point will be a breakpoint. Some 
thresholds are computed: LNTSH1 is such that at any point on the boundary, 
if the distance to the next point is greater than LNTSH1, then this seg- 
ment is automatically taken to be a straight line piece. If the said dis- 
tnace is less than LNTSH2 then it is assumed that this will not be a 
straight line piece, but rather part of some curve piece. The other cases 
are resolved by some additional formulae: consider point i with associ- 



ated distance d. and inside angle a! (outside angle a. = it - a!), 
l i l l 
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(x.y.) 



(x i+iyi+i ) 




We will consider the segment (x ± y i ) - ^n-i^n-i^ to be a straight line 



piece if d. > LNTSH1 or if 
1 



[d. > LNTSH2 /\ (a.,- " d. > CRVF v a? • d. > CRVB)]. 
l l+l l 11 

Any curved piece is extended until an inflection point or straight line 
is reached. From the above formula, it is clear that straight lines are 
usually found by the local property of "sharp corner". The formula involves 
the two endpoints of a straight line in slightly different ways because it 
was found to be efficacious in discerning certain odd types of corners while 
using a one-way scan of the points from ANGLES. For example 




if 



or 




Experiments were carried out to determine suitable values for the above 
thresholds, with the following results: 
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LNTSH1 = 0.25 * DIAMETER 
LNTSH2 = 0.03 * DIAMETER 
CRVB = 0.022 * DIAMETER 
CRVF = 0.02068 * DIAMETER 

The curved pieces output by SLPCS may be represented as some segment 
of a conic (general second degree curve) since there are no inflection 
points in such curves (some pathological cases do arise however, such as 
egg — shaped closed curves and spiral sections, but for the moment these 
are ignored, and probably will not be parts of viewed objects anyway). 
BSTCRV accepts any piece from the output of SLPCS, including straight 
lines which it leaves unaltered, and returns a list of eight members repre- 
senting a conic-curve segment. That is, for straight lines 



(0 (x lYl ) (x 2 y 2 )) BSTCRV > (0 (x lYl ) (x 2 y 2 )) 



and for curves which can be approximated by .a circular or, elliptical arc 

(IND (x lYl ) ... (x 2 y 2 )) BSTCRV, (IND ! ELLPSE f P £\ A a P'lB) 

where 

IND is 1 or 2 for curved out or in 

T ELLPSE ? is the type of curve, as opposed to hyberbolic 

P(3 is a point representing the center of the ellipse 

Xl is a LISP number for the length of the principal 

semi-axis 

X 2 is for the secondary semi-axis 
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a is an angle -tt/2 < a -^ tt/2 measuring the tilt of the 

principal axis from the horizontal, about the center 
of the ellipse. 

Pi, 6 specify the segment of the curve in the following 

manner: "beginning at point Pi and using Pq as cen- 
ter point, take B radians about Pq, in the clock- 
wise direction if IND = 1 otherwise in the counter- 
clockwise direction. This is the relevant segment. 

Clearly for a circle li = ^2, and for any complete circle or ellipse, 

3 = 2tt and Pi is any point on the curve. 

Example 




For second-degree curves which are not ellipses, i.e. hyperbolae or 
parabolae we simply return the general-equation coefficients and the 
first and last point on the arc. 

(IND (xiyx) ... (x y 5) BSTCRV^ (IND 'HYPER' (A B) C D E (x^) (x y )) 
where the best fit curve is 

Ax 2 + Bxy + Cy 2 +Dx + Ey = 1 



rpippjiyww 
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Objects which have no inflection points 
returned by SLPCKS as all erne piece. A* 
objects fall into this category, b»* they 
Ixcept for a few immiu, BSTdV can be 
by MM* of Naf€*S — 8XL being the 
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, egg-shaped 
-fit to a conic, 
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Appendix I 

Atoms Used as Free Variables 

Throughout the Vision Programs 



PI 
TWOPI 



Set to 7T and 2tt by a top level call to SETQ, these 
atoms provide convenient reference to some familiar 
numerical values. 



CTOL 
ETOL 



BSTCRV requires the mean -square error in an approximation 
to be less than some tolerance. Since circles are intrin- 
sically more simple than general conies, a circular fit is 
first examined: if the computed circle has center (a,b) 
and radius v, then it is accepted iff 



i|^ 2 + fei^) 2 _ (1+ ^^)" 

n 4^ 1 1 r2 r z r z 



< CTOL 



If the circle is not accepted, a general second-degree 
curve fit is tried and accepted iff 

2 



u 



i=l 



Ax 



? + Bx.y. + Cy? + Dx. + Ey . -1 
1 i/i ; i i ; i 



< ETOL 



See Appendix II for the formulae by which the best fit 

curve coefficients are determined. 

Suitable values for these tolerances are ETOL = .01 

CTOL = .0015 
but some applications may want to vary them. 



mmmmmmm mmmmm»^mi>»mm>^mmm WMUJIiifPHUu wmwwmvvvvpwm 



r*s 
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POLYStS determines a straight lis* approximation to a 
curve by initially pUfctof »P WL p#li*a aad extending 
by H2 points at a tine until the normal ixed mean-square 
error of the beet-fit etraaght line to these points is 
greater than IK.. "IMf-flf" meant leaat-aquarea beet 
fit in which 

— ■ > (ax. + by. -e} 2 ia minimised by appropriate choice 

f. i-1 X * 

of a,b, and c. 

More experimentation should be done te find suitable 
values for these parameters — a good lirst guess might 
be X9L - .05 

m. - 4 

*2 - 4 
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Appendix II 

BEST FIT CIRCLE 

For a data set ((x^) ... (x y )] the best fit circle, with equa- 

9 9 

tion ax' 1 + bx + ay* 1 + cy = 1.0, is the one which minimizes 

.2 



i i 



ax: + bx. + ay? + cy . - 1.0 



n 1=1 

Differentiating the error expression above with respect to a, b, and c, 
leads to three simultaneous linear equations, which then can easily be 
solved. The center of the circle is at 



(- -Tr-f - tt-) and it has radius x / — + A 2 + B 2 
za la. V ,a 



BEST FIT CONIC 

As above, we solve five simultanious linear equations in order to 
minimize 



-£[ax 2 + bxy + cy 2 + dx + ey -1.0] 2 . 



Actually, a parametric form of the second-degree equation would be more 

revealing, 

ax 2 + bxy + cy + dx + ey = f . 

One must place some additional constraint, however, since there ^are only 

five degrees of freedom in a conic curve, but the parametric equation has 

six unknowns. Imposing 

\la 2 + b 2 + c 2 + d 2 + e 2 + f 2 = 1 
looks plausible, but the resulting simultaneous equations are horribly 
non-linear, so the much simpler condition f = 1 has been chosen. Cases in 
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which this latter assumption leads to numerical nonsense rarely occur 
in practice. 

DETERMINATION OF ELLIPSE PARAMETERS 

If b - 4ac < 0, then the curve represented by 

ax z + bxy + cy 2 + dx + ey = 1.0 
is an ellipse. Otherwise it is a parabola or hyperbola. If the whole 
plane is rotated about the origin through an angle of -a radians, the 
axes of the ellipse will be parallel to the x- and y-axes. Thus consider 
the curve 

Ax 2 + Cy 2 + Dx + Ey = 1.0 

where 

1 -1 b f 1 if x > 

a = sign (b) • sign(a-c) • — tan | 1 , sign(x) = < 

1 a ° (^l if x < 

A = a cos a + b cosa sina + c sin a 

B = b(cos a - sin a) + 2(c-a) sina cosa = 

C = a'sin a - b-cosa • sina + c cos a 

D = d cosa + e sina 

E = -d sina + e cosa 

D E 
This represents an ellipse with center at (- — - y, - ^r) , whose horizontal 

semi-axis length is 

/ - \/4AC 2 + D 2 C 2 + ACE 2 " 
2AC 

and whose vertical semi-axis length is^/2 = 77 '^1- 









Thus the origin h*a 






( \ [-8| cam + ^rnimm h § C--j^ 



Reference: Thoaae, G.B. Jr* 
Wesley U940). 
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